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Abstract 



We calculate the relic density (O^o/i^) of neutralino dark matter in su- 
pergravity models with radiative electroweak symmetry breaking, for values 
of the Higgs sector parameter tan/3 = V2/V1 between the two infrared fixed 
points of the top quark Yukawa coupling, tan/3 ~ 1.8 and tan/3 ~ 56. For 
tan/3 < 10, the CP-odd Higgs pseudoscalar {A^) and the heavier CP-even 
Higgs scalar (H^) are much heavier than the lightest neutralino (x?); and 
the annihilation cross section of Xi pairs is significantly enhanced only in 
the vicinity of the lighter CP-even Higgs boson (h^) and the Z boson poles. 
For tan/3 ^ 40, niA and mn become comparable to 2m^o, and the neu- 
tralino annihilation cross section can be significantly enhanced by the poles 
of the and the as well. For a cosmologically interesting relic density, 
0.1 ^ Q^oh^ ^ 0.5, we find that the supergravity parameter space is tightly 
constrained. The mass spectra for supersymmetric particles are presented for 
unification mass parameters in the cosmologically interesting region. 



Typeset using REVT^ 



I. INTRODUCTION 



A super symmetry (SUSY) between fermions and bosons provides a natural explanation 
of the Higgs electroweak symmetry breaking (EWSB) mechanism in the framework of a 
grand unified theory (GUT). The evolution of renormalization group equations (RGEs) |l| 
with the particle content of the minimal supersymmetric standard model (MSSM) is 
consistent with a grand unified scale at Mgut ~ 2 x 10^^ GeV and an effective SUSY mass 
scale in the range Mz < Msusy ~ 1 TeV. With a large top quark Yukawa coupling to a 
Higgs boson at the GUT scale, needed for h—r Yukawa unification radiative corrections 
can drive the corresponding Higgs boson mass squared parameter negative, spontaneously 
breaking the electroweak symmetry and naturally explaining the origin of the electroweak 
scale. In the minimal supersymmetric GUT with a large top quark Yukawa coupling (y^), 
there is an infrared fixed point (IRFP) at low tan/?, and the top quark mass is cor- 

respondingly predicted to be rrit = 200 GeV sin /5 0. At high tan/3, another quasi-IRFP 
solution {dYt/dt ~ 0) exists at tan/? ~ 60 For mt = 175 GeV, the two IRFPs of top 

quark Yukawa coupling appear at tan/5 ~ 1.8 and tan/? ~ 56 |Q. 

Supersymmetry extends the Poincare symmetry to new dimensions of space and time 
and local gauge invariance with SUSY includes gravity. In supergravity (SUGRA) models 
P], it is assumed that supersymmetry is broken in a hidden sector with SUSY breaking 
communicated to the observable sector through gravitational interactions, leading naturally 
but not necessarily |jlOl to a common scalar mass (mg), a common gaugino mass {mi/2), 
a common trilinear coupling {Aq) and a common bilinear coupling (Bq) at the GUT scale. 
Through minimization of the Higgs potential, the B parameter and magnitude of the super- 
potential Higgs mixing parameter /i are related to the ratio of vacuum expectation values 
(VEVs) of Higgs fields (tan/5 = V2/V1), and the mass of the Z boson (Mz). The SUSY 
particle masses and couplings at the weak scale can be predicted by the evolution of RGEs 
from the unification scale 



In our analysis, we evaluate SUSY mass spectra and couplings in the minimal supergrav- 
ity model with four parameters: mo, mi/2, Aq and tan/?; and the sign of the Higgs mixing 
parameter fi. The neutralino relic density is only slightly affected by a change in the Aq 
parameter since Aq mainly affects the masses of third generation sfermions. Therefore, for 
simplicity we take Aq = in most of our calculations. 

The mass matrix of the neutralinos in the weak eigenstates {B, W3, Hi, H2) has the 
following form |^ 

I Ml -M2Cos/3sin6'iy Mz sin /5 sin 6^^^ \ 

Mn = 



V 



" yy ---Zy Vy 

M2 Mz COS /? COS 6*^^ — Mz sin/5cos6'^y 

-Mz cos (3 sin Mz cos [3 cos Q n 

Mz sin P sin 9yy —Mz sin /5 cos 9yy /i J 



(1) 



This matrix is symmetric and can be diagonalized by a single matrix to obtain the mass 
eigenvalues and eigenvectors. The form of Eq. (|l]) establishes our sign convention for /i. 
Recent measurements of the b s'j decay rate by the CLEO [|T2[ and LEP collaborations 



p!3| place constraints on the parameter space of the minimal supergravity model |jT^. It was 



found that 6 — > 57 excludes most of the minimal supergravity (mSUGRA) parameter space 
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when tan/3 is large and /i > [jl^. In this article, we present results for both /i > and 
/i < 0. 

We calculate the masses and couplings in the Higgs sector with one loop corrections from 
both the top and the bottom Yukawa interactions in the RGE-improved one-loop effective 
potential at the scale Q = ^Jfoijni^. With this scale choice, the numerical value of 
the CP-odd Higgs boson mass (m^) at large tan/5 is relatively insensitive to the exact 
scale and the loop corrections to Ma are small compared to the tree level contribution. In 
addition, when this high scale is used, the RGE improved one-loop corrections approximately 
reproduce the dominant two loop perturvative calculation of the mass of the lighter CP-even 
Higgs scalar ivrih). Our numerical values of vrth are very close to the results of Ref. where 
somewhat different scales higher than Mz have been adopted in evaluating the effective 
potential. 

The matter density of the Universe p is described in terms of a relative density Vt = p/ pc 
with pc = 3Hq/87tGn — 1-88 x 10~^^/i^ g/cm^ the critical density to close the Universe. 
Here, Hq is the Hubble constant, h = ifo/(100 km sec~^ Mpc~^), and Gn is Newton's 
gravitational constant. 

There is compelling evidence from astronomical observations for the existence of dark 
matter |TB[ . Studies on clusters of galaxies and large scale structure suggest that the matter 
density in the Universe should be at least 20% of the critical density {Qm ~ 0.2) |jl9[, but 
the big-bang nucleosynthesis and the measured primordial abundance of helium, deuterium 
and lithium constrain the baryonic density to 0.01 ^ fi^/i^ ^ 0.03 ||20|. The anisotropy in 



the cosmic microwave background radiation measured by the Cosmic Background Explorer 
(COBE) suggests that at least 60% of the dark matter should be cold (nonrelativistic) 
P 1|| . In supersymmetric theories with a conserved i?-parity[|, the lightest supersymmetric 
particle (LSP) cannot decay into normal particles and the LSP is an attractive candidate for 
cosmological dark matter []22| , |23| . In most of the supergravity parameter space, the lightest 
neutrahno (xi) is the LSP 

Inflationary models usually require f2 = 1 for the Universe p2|, although models with 



i7 < 1 have recently been proposed P4 
converging to ^ 0.6 — 0.7 p 



Recent measurements on the Hubble constant are 



Therefore, we conservatively consider the cosmologically 



interesting region for f2^o to be 



0.1 < n^oh^ < 0.5. 



(2) 



Following the classic papers of Zel'dovich, Chiu, and Lee and Weinberg PU[, many studies 



of the neutralino relic density in supergravity models have been made with universal [p^^ -|30 



and nonuniversal |31| soft breaking masses. Recent calculations have taken into account 
threshold effects and integration over Breit-Wigner poles [^,0. Previous analyses have 
focused on the small to moderate tan/3 range, tan/3 ^ 25, but large tan/3 is equally of 



interest. Some theoretical models of fermion masses require a large value tan/3 |^ 



In this article, we make more complete calculations of the neutralino relic density over the 
full range of tan /3 and SUGRA parameter space with universal GUT scale boundary condi- 
tions. We employ the calculational methods of Ref. 130] , with some significant improvements. 



R = +1 for the SM particles and Higgs bosons and R = —1 for their superpartners. 
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The neutralino annihilation cross section is evaluated with helicity amplitude techniques to 
eliminate any uncertainty from a neutralino velocity expansion. The neutralino relic density 
is calculated with relativistic Boltzmann averaging, neutralino annihilation threshold effects 
and Breit-Wigner poles. We introduce a transformation to improve the numerical integra- 
tion over Breit-Wigner poles and use appropriate extension of the RGEs for large tan/? 
solutions. We determine the constraints on mi/2 and tan/3 implied by Eq. (||). We also 
impose the constraints from the chargino search at LEP 2 [^. The effects of the common 
trilinear coupling (Aq) on At, and the neutralino relic density Q^oh'^ are also studied. 



II. RELIC DENSITY OF THE LIGHTEST NEUTRALINO 

The present mass density of the lightest neutralino is inversely proportional to the an- 
nihilation cross section of Xi pairs. When kinematically allowed, the neutralino pairs an- 
nihilate into fermion pairs (//), gauge boson pairs and ZZ), Higgs boson pairs[| 
H°H°, A°A°, h^H^ h^A^, H°A° and H+R-) and associated pairs of gauge and Higgs 
bosons {Zh^, ZH^, ZA^, and W^H^) through s, t, and u channel diagrams. The most con- 
venient approach to calculate the annihilation cross section including interferences among 
all diagrams is to calculate the amplitude of each diagram with the helicity amplitude for- 
malism]^ then numerically evaluate the matrix element squared with a computer. Therefore, 
following Ref . , we employ the HELAS package |]36| to calculate the helicity amplitudes. 



The annihilation cross section of Xi pairs depends on its mass {m^o) which we find can be 
empirically expressed in GeV units as a function of the GUT scale gaugino mass mi/2 and 
sin 2/3 as 

/i > : m^o ~ 0.448mi/2 + 11.7 sin 2/3 - 10.4, 

/i < : m^o ~ 0.452mi/2 + 4.68 sin 2(3 - 12.9, (3) 

for 100 GeV ^ ~ 1000 GeV and all tan/3 for which perturbative RGE solutions exist. 
These neutralino mass formulas hold to an accuracy of ^ 3%. 

In the regions of Breit-Wigner resonance poles the near-singular amplitudes make the 
integration unstable. To greatly improve the efficiency and accuracy in calculating the an- 
nihilation cross section for processes with a Breit-Wigner resonance, we introduce a trans- 
formation s — = Mr tan 6 such that 1571 



ds f d9 

(4) 



(s - M2)2 + (Mr)2 J MT 

and scale 9 to the range (0,1). 

The time evolution of the number density n{t) of weakly interacting mass particles is 
described by the Boltzmann equation p6| 



2There are five Higgs bosons: two CP-even (lighter) and (heavier), one CP-odd A^, and a 
pair of singly charged H^. 



The helicity amplitudes of all relevant processes in the non-relativistic limit are given in Ref. |28|. 
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dt 



{av) [n^ - n|] 



(5) 



where H = IMgl^^T^/Mpi is the Hubble expansion rate with (yf^, ~ 81 the effective number 
of relativistic degrees of freedom, Mpi = 1.22 x 10^^ is the Planck mass, (av) is the thermally 
averaged cross section times velocity, v is the relative velocity and a is the annihilation cross 
section, and is the number density at thermal equilibrium. 

In the early Universe, when the temperature T ^ "^x?' ^iSP existed abundantly 
in thermal equilibrium with the LSP annihilations into lighter particles balanced by pair 
production. Deviation from the thermal equilibrium began when the temperature reached 
the freeze-out temperature (Tf ~ m^o/20). After the temperature dropped well below Tf, 
the annihilation rate became equal to the expansion rate and n^o = H/ (crv). The resulting 
relic density can be expressed in terms of the critical density by p2 



m 



(6) 



where n^o is the number density of Xi- 

The thermally averaged cross section times velocity is 



{av) (T) 



Je-Ei/Te-E,/T^3p^^3p^ 



(7) 



where pi (Ei) and p2 (-E'2) are the momentum and energy of the two colliding particles in 
the cosmic, co-moving frame of reference, and T is the temperature. This expression has 



been simplified to a one-dimensional integral |33] of the form 



(av) (x) 



1 f°° 

727T\ / daa{a)a'^{a'^ — 4:)Ki 

^9 ( Z ) -^2 



(8) 



where x = a = ^/s is the subprocess energy, and the Ki are modified Bessel 

functions of order i. 

To calculate the neutralino relic density, it is necessary to first determine the scaled 
freeze-out temperature xp = ^xll^^f- '^^^ standard procedure is to iteratively solve the 
freeze out relation 



xJ = loe 



m^o / 45 1/2 



(9) 



starting with xp = The relic density at the present temperature (Tq) is then obtained 
from 



Pin) 



8.0992 X 10-47 GeV^' 



(10) 



where 



Pin) ^ 1.66 X ^ T^y^ 



Mpi \ n I Jq'^ {av)dx 



(11) 
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The ratio T^o/T^ is the reheating factor of the neutrahno temperature (T^o) compared 
to the microwave background temperature (T^). The integration in Eq. ([TTl) is evaluated 
with the modified Bessel functions expanded as power series in x, and then integrated over 
X. The result is 



Xp 



(cTf ) dx 



daa{a)a^/^{a^ - 4)F(a) 



(12) 



where 30 



F{a) 




e 



Xp — 2{a — 2)^/xF 



(13) 



Essentially, all of the contribution to the integral in Eq. (|T2|) comes from a < 2.5. 

To generate the SUSY particle mass spectrum and couplings for the neutrahno relic 
density calculation, we run the gauge coupling RGEs from the weak scale up to a high 
energy scaleQ to fix the unification scale Mgut, and the unified gauge coupling oqut- Then 
we run all the RGEs from Mqut to the weak scale of Mz- The SUSY mass scale (Msusy) 
is chosen to be 2 TeV; closely similar results are obtained with a 1 TeV mass scale. 

In Figure 1, we present masses, in the caseQ /i > 0, for the lightest neutrahno (xi)) 
the lighter top squark (ti), the lighter tau slepton (fi), and two neutral Higgs bosons: the 
lighter CP-even {h^) and the CP-odd {A^). To a good approximation, the mass of the 
lightest chargino {xt) is about twice m^o and the mass of the heavier CP-even Higgs boson 
(i7°) is close to uia [1,0] • For rriQ ~ 100 GeV and mi/2 ^ 400 GeV, the mass of fi can 
become smaller than m^o so such regions are theoretically excluded. Also shown in Fig. 1 are 
the regions that do not satisfy the following theoretical requirements: electroweak symmetry 
breaking (EWSB), tachyon free, and the lightest neutrahno as the LSP. The region excluded 
by the > 85 GeV limit from the chargino search |^ at LEP 2 is indicated. Masses for 
yU < are shown in Fig. 2. 

There are several interesting aspects to note in Figs. 1 and 2: 
(i) For tan/3 ~ 1.8, n > generates heavier neutralinos and charginos but lighter than 



^ For simplicity, we employed the one-loop RGEs in our analysis since the two-loop modifications 
to the SUSY mass spectra are only a few percent of the one-loop RGE results |^ . We did not include 
heavy particle threshold effects in the RGE evolution, which may introduce small corrections when 
tan/? is large [pS|] . 

^ Our convention for the sign of fi follows that of Ref. and is opposite to that of Ref. see 
Eq. dlD. 
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/i < for the same mi/2 and mo. 

(ii) An increase in tan/3 leads to a larger rrih but a reduction in m^o, Tn^±, ttla and ttlh- 

(iii) For tan/3 ~ 50, m^i ~ m^^ can become comparable to fn^o. 

(iv) Increasing mo raises m^, "nT.^ and the masses of other scalars significantly. 

The relic density of the neutralino dark matter {fl^oh'^) for the case /x > is presented 
in Fig. 3 versus mi/2 for several values of tuq and tan/3. Figure 4 shows similar contours for 
/i < 0. When 2m^o is close to the mass of Z, h^, H^, or A^, the cross section is significantly 
enhanced by the Breit-Wigner resonance pole, and a corresponding dip appears in the relic 
density distribution. For fi > 0, tan/? = 1.8 and mo = 100 GeV, there are two dips, the first 
one occurs at the pole (mi/2 = 55 GeV, m^o = 26 GeV and rrih = 52 GeV), and the second 
one appears at the Z pole (mi/2 — 100 GeV and m^o = 45 GeV ~ |M^). An increase in mo 
raises the mass of Higgs bosons but does not affect the mass of Xi- Therefore, comparing 
the mo = 100 GeV results to those of mo = 500 GeV, the Z-pole dips occur at about the 
same mi/2 as above, but the h-pole dips appear at larger mi/2 and m^o. Increasing tan/3 
slightly raises rrih, slightly lowers m^o, and greatly reduces m^ and mn- For tan/3 ~ 45, 
extra dips appear in Figs. 3(c) and 4(c). These resonance dips are associated with the 
and masses and occur where 2m^o ~ m^ ~ rriH- For tan/3 < 10 and mo > 200 GeV, 
the annihilation cross section is too small for a cosmo logically interesting fl^oh'^, while for 
tan /3 ~ 45, mo must be close to 1000 GeV for mi/2 ^ 750 GeV and be larger than 400 GeV 
for mi/2 ^ 750 GeV to obtain an interesting relic density. For tan/3 > 50, m^ < ^m^o, and 
there are no A^ and resonant contributions; an acceptable flh'^ is found for mo ^ 400 
GeV and 400 GeV < mi/2 < 800 GeV. 



III. CONSTRAINTS ON THE SUGRA PARAMETER SPACE 

To show the constraints of Eq. (H) on the SUGRA parameter space, we present contours 
of f^x?^^ ~ (mi/2,tan/3) plane in Fig. 5 for /i > 0. Also shown are the 

regions that do not satisfy the theoretical requirements (electroweak symmetry breaking, 
tachyon free, and lightest neutralino as the LSP) and the region excluded by the chargino 



search (m + < 85 GeV) at LEP 2 [^. Figure 6 shows similar contours for n < 0. 



We summarize below the central features of these results. 
If mo is close to 100 GeV, 

• Most of the (mi/2,tan/3) plane with tan/3 ^ 40 is excluded by the above mentioned 
theoretical requirements. 

• The chargino search at LEP 2 excludes the region where mi/2 ^ 100 GeV for n > 

and mi/2 < 120 GeV for jj, <0. 

• The cosmologically interesting region is 100 GeV ^ mi/2 ~ 400 GeV and tan/3 ^ 25 
for either sign of fi. 

If mo is close to 500 GeV, 

• Most of the (mi/2,tan/3) plane is theoretically acceptable for tan/3 ^ 50 and mi/2 ^ 
120 GeV. 
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• The LEP 2 chargino search excludes (i) mi/2 ~ 140 GeV for tan/9 ^ 10, and (ii) 
mi/2 ^ 80 GeV (/x > 0) or mi/2 ^ 120 GeV {pi < 0) for tan/? - 1.8. 

• The cosmologically interesting regions he in two narrow bands. 

Contours of fl^oh'^ = 0.1 and 0.5 in the (mi/2, mo) plane are presented in Fig. 7 for /i, > 
with tan/3 = 1.8 and 50. Also shown are the theoretically excluded regions. Figure 8 shows 
similar contours for /i < 0. Salient features of these results are summarized in the following. 
If tan/3 is close to 1.8, 

• Most of the (mi/2,mo) parameter space is theoretically acceptable. 

• The chargino search at LEP 2 excludes the region where mi/2 ^ 80 GeV for > 
and mi/2 ^110 GeV for < 0. 

• Most of the cosmologically interesting region is 80 GeV ^ mi/2 ^ 450 GeV and 
mo < 200 GeV. 

If tan P is close to 50, 

• The theoretically acceptable region in the (mi/2,mo) plane is constrained to have 
mo > 160 GeV and mi/2 > 150 GeV for tan /3 ~ 50. 

• The LEP 2 chargino search excludes (i) the region with mi/2 ^ 125 GeV for > or 
(ii) the region with mi/2 ^ 135 GeV for < 0, which is already inside the theoretically 
excluded region. 

• The cosmologically interesting region lies in a band with (i) 475 GeV ^ mi/2 ^ 800 
GeV for // > 0, or (ii) 500 GeV < mi/2 < 840 GeV for // < 0, and mo > 300 GeV. 

To correlate the acceptable values of neutralino relic density to the masses of SUSY 
particles at the weak scale, we present ft^oh^ and representative SUSY mass spectra versus 
mi/2 in Fig. 9, for > with (a) tan/3 = 1.8, mo = 150 GeV and (b) tan/3 = 50, mo = 600 
GeV. Figure 10 shows similar contours for fi < 0. Also shown are the theoretically excluded 
regions and the excluded region from the chargino search at LEP 2. 

The neutralino rehc density and representative SUSY mass spectrum versus mo are 
presented in Fig. 11, for > with (a) tan/3 = 1.8, mi/2 = 150 GeV and (b) tan/3 = 
50, mi/2 = 500 GeV. Figure 12 shows similar contours for fi < 0. Also shown are the 
theoretically excluded regions. 

The foregoing results were based on the GUT scale trilinear coupling choice Aq = 0. To 
examine the sensitivity of the results to ^40, we present the neutralino relic density versus 
Aq in Fig. 13 for various values of mo = mi/2 and tan/3. For tan/3 < 10, Q^o/i^ is almost 
independent of Aq. For tan j3 ~ 40, the relic density is reduced by a positive Aq but enhanced 
by a negative Aq. Significant impact of Aq on occurs only when both mo and mi/2 

are small and tan,/3 is large. For mo = mi/2 = 200 GeV and tan/3 = 40 the theoretically 
allowed range of Aq is -380 GeV < Aq < 590 GeV. The value of fi^o/i^ with Aq = 590 
GeV (-380 GeV) is about 0.22 (2.4) times the value with Aq = 0. For mo = mi/2 = 800 
GeV and tan/3 = 40 the value of fi^o/i^ with Aq = 1000 GeV (-1000 GeV) is about 0.58 
(1.7) times the value with Aq — 0. 
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The dependence of ^^l^^ occurs because the masses of top squarks are sensitive 

to At, and the t-channel diagrams involving the stops make important contribution to neu- 
trahno annihilation cross section. In figures 14 (/i > 0) and 15 (/i < 0), we show the mass 
of the lighter top squark (rn^-J and the trilinear coupling At versus Aq for various values 
of mo = ?Tii/2 and tan/5. The trilinear coupling At is always negative and it is nearly a 
linear function of Aq. For tan/5 close to the top Yukawa infrared fixed point, tan/5 ~ 1.8, 
At approaches an infrared fixed point vaule as well. At ~ — 1.7mi/2, to about 30% accuracy. 



independent of the value of Aq [40 



IV. CONCLUSIONS 

The existence of dark matter in the Universe provides a potentially important bridge 
between particle physics and cosmology. Within SUGRA GUTs, the lightest neutralino 
is the most likely dark matter candidate. We evaluated the neutralino relic density using 
helicity amplitude methods including all relevant diagrams. We applied a transformation 
to improve the efficiency and accuracy of annihilation cross section calculations in regions 
with Breit-Wigner resonance poles. Requiring that the neutralino relic density should be in 
the cosmologically interesting region, we were able to place tight constraints on the SUGRA 
parameter space, especially in the plane of mi/2 versus tan/5, since the mass of the lightest 
neutralino depends mainly on these two parameters. We derived empirical formulas relating 
m^o to mi/2 and sin/5. The cosmologically interesting regions of the parameter space with 
tan /5 close to the top Yukawa infrared fixed points are found to be 

tan/5 = 1.8 : 80 GeV < mi/2 < 450 GeV and mo < 200 GeV 
tan/5 = 50 : 500 GeV < mi/2 < 800 GeV and mo > 300 GeV 

where the high tan /5 result is based on Aq = 0. Both regions are nearly independent of the 
sign of /i. We presented expectations for the SUSY particle mass spectra corresponding to 
these regions of parameter space. For the IRFP at tan/5 = 1.8, the lighter scalar tau (fi) 
and the lightest scalar neutrino {u) can be relatively light 200 GeV) in the cosmologically 
interesting region. 

For tan/5 ^ 10, the neutralino relic density Q^oh'^ is almost independent of the trilinear 
couplings Aq. For tan/5 ~ 40, the rehc density is reduced by a positive Aq while enhanced 
by a negative Aq. The value of Aq significantly affects ^x?^^ only when tan/5 is large and 
both mo and mi/2 are small. 

Supernovae ETI] , cosmic microwave background P2| and other astronomical observations 



will soon provide precise information on the Hubble parameter h and the cold dark matter 
component of Q. This information will pinpoint the SUGRA GUT parameters, which in 
turn predict the SUSY particle masses in these models. 
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FIG. 1. Masses of Xii ti at the mass scale of Mz and masses of h^, and at the mass 
scale Q = ^m^-^ m^^ , versus mi/2, with MsusY = 2 TeV and > for (a) tan/3 = 1.8, mo = 100 
GeV, (b) tan/3 = 1.8, mo = 500 GeV, (c) tan /? = 50, mo = 300 GeV, and (d) tan /? = 50, mo = 600 
GeV. The shaded regions are excluded by theoretical requirements (EWSB, tachyons, LSP), or the 
chargino search at LEP 2. 
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FIG. 2. The same as in Fig. 1, except that < 0. 
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FIG. 3. ^x?^^ versus mi/2 ^or fi > 0, with representative values of mo and (a) tan/3 = 1.8, 
(b) tan/3 = 10, (c) tan/3 = 45, and (d) tan/3 = 50. The shaded regions denote the parts of 
the parameter space (i) producing ^-^oh'^ < 0.1 or ^^^o/i^ > 0.5, (ii) excluded by theoretical 
requirements, or (iii) excluded by the chargino search at LEP 2. 
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FIG. 4. The same as in Fig. 3, except that fi < 0. 
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FIG. 5. Contours of ^^oh'^ = 0.1 and 0.5 in the (mi/2)tan/3) plane, for /i > 0, (a) mo = 100 
GeV and (b) niQ = 500 GeV. The shaded regions denote the parts of the parameter space (i) pro- 
ducing ri^o/i^ < 0.1 or ^^oh'^ > 0.5, (ii) excluded by theoretical requirements, or (iii) excluded by 
the chargino search at LEP 2. 
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FIG. 6. The same as in Fig. 5, except that /U < 0. 
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and (b) tan f3 = 50. The shaded regions denote the parts of the parameter space (i) producing 
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FIG. 9. The neutrahno rehc density and representative SUSY mass spectrum versus mi/2 
for ^ > with (a) tan/3 = 1.8, mo = 150 GeV and (b) tan/3 = 50, mo = 600 GeV. In Figs. 
9-12, v is the hghtest scalar neutrino. The shaded regions denote the parts of the parameter 
space (i) producing ^^oh? < 0.1 or ^^oh"^ > 0.5, (ii) excluded by theoretical requirements, or 
(iii) excluded by the chargino search at LEP 2. 
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FIG. 10. The same as in Fig. 9, except that /x < 0. 
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(a) tan/S=1.8, mi/2=150 GeV, ix>0 (b) tan|8=50, mi/g=500 GeV, //>0 
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FIG. 11. The neutralino relic density and representative SUSY mass spectrum versus txiq for 
^ > with (a) tan/3 = 1.8, mi/2 = 150 GeV and (b) tan/? = 50, m;^/2 = 500 GeV. The shaded 
regions denote the parts of the parameter space (i) producing ^^oh'^ < 0.1 or ^^o^h"^ > 0.5, or 
(ii) excluded by theoretical requirements. 
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(a) tan/S=1.8, mi/g=150 GeV, /i,<0 (b) tan/S=50, mi/2=500 GeV, /j.<0 
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FIG. 12. The same as in Fig. 11, except that ^ < 0. 
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FIG. 13. The neutralino relic density versus the GUT scalar trilinear coupling Aq for (a) // > 
and tan/3 = 1.8, (b) /i > and tan/3 = 40, (c) /x < and tan/3 = 1.8, (d) /U < and tan/3 = 40, 
and various values of mo = mi/2- 
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FIG. 14. The mass of the hghter top squark versus Aq for (a) tan/? = 1.8, (b) tan/3 = 40, and 
various values of mo = mi/2 and n > 0. Also shown is the soft breaking trilinear coupling At at 
the weak scale versus the GUT scale input Aq for (c) tan/3 = 1.8, (d) tan/3 = 40, 
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FIG. 15. The same as in Fig. 14, except that /x < 0. 
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